Condensate fraction in metallic 
superconductors and ultracold atomic vapors 

Luca Salasnich 

CNR and CNISM, Dipartimento di Fisica "Galileo Galilei" , 
Universita di Padova, Via Marzolo 8, 35131 Padova, Italy 
e-mail: luca.salasnich@cnr.it 



Abstract 

We investigate the condensate density and the condensate fraction 
of conduction electrons in weak-coupling superconductors by using the 
BCS theory and the concept of off-diagonal-long-range-order. We dis- 
cuss the analytical formula of the zero-temperature condensate density 
of Cooper pairs as a function of Debye frequency and energy gap, and 
calculate the condensate fraction for some metals. We study the den- 
sity of Cooper pairs also at finite temperature showing its connection 
with the gap order parameter and the effects of the electron-phonon 
coupling. Finally, we analyze similarities and differences between su- 
perconductors and ultracold Fermi atoms in the determination of their 
condensate density by using the BCS theory. 

PACS numbers: 74.20.Fg; 74.70.Aq; 03.75. Ss. 



1 Introduction 

The condensate fraction of fermionic alkali-metal atoms has been recently 
investigated [D [21 S i] by using extended BCS (EBCS) equations [5l il [3 El 
m [To] from the BCS regime of Cooper-pairs to the BEC regime of molecular 
dimers [H [21 [3]. In particular, we have found [1] a remarkable agreement 
between this simple mean-field theory and the experimental results [HI IT2] . 
These results indicate the presence of a relevant fraction of condensed pairs of 
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^Li atoms also on the BCS side of the Feshbach resonance. Monte Carlo cal- 
culations [I3] have shown that the zero-temperature mean-field predictions 
[U [2] slightly overestimate the condensed fraction of Fermi pairs. Very re- 
cently it has been reported |T3] an accurate measurement of the temperature 
dependence of the condensate fraction for a fermion pair condensate of ^Li 
atoms near the unitarity limit of the BCS-BEC crossover. Also these new ex- 
perimental data [11] are in agreement with mean-field theoretical predictions 
at finite temperature [3]. 

In superfiuids made of ultracold atoms, the inter-atomic interaction is 
attractive for all fermions of the system [6] . On the contrary, in metallic su- 
perconductors there is an attractive interaction between fermions only near 
the Fermi surface [TSl [IS]. As a consequence, the condensate fraction of 
metallic superconductors has distinctive properties with respect to those of 
atomic superfiuids. Despite the BCS theory is 52 years old \1T\, the conden- 
sate fraction of Cooper pairs in superconductors has been considered only in 
few papers |lHl EOl [191 12l] and in the recent book of Leggett |15]. In fact, 
in superconductors the condensate fraction has never been measured: only 
very recently Chakravarty and Kee have proposed to measure it by using 
magnetic neutron scattering [21] . 

In this paper we analyze in detail the condensate density of conduction 
electrons in weak-coupling superconductors at zero and finite temperature 
by using BCS theory [T7] and the concept of off- diagonal-long-range-order 
[TH[ [22] • For the first time, we calculate explicitly the density of electronic 
Cooper pairs and the condensate fraction for various metals and show its 
dependence on the Debye frequency, the electron-phonon interaction and the 
energy gap. Another novelty of this paper is the analytical and numerical 
investigation of the temperature dependence of the condensate fraction, for 
which we find a power-law behavior. Finally, we compare of the BEC theory 
of superconductors with the extended BEC theory of ultracold Fermi atoms 
for obtaining the condensate density and the condensate fraction. 



2 BCS theory and ODLRO 

The BCS Lagrangian density of conduction electrons with spin cr = f, J, near 
the Fermi surface is given by 

^ = J2^t f - e( V) + ^^^. + 9 ^t^v^;+v^4^t , (1) 
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where ip„{r,t) is the electronic field operator which satisfies the familiar 
equal-time anti-commutation rules of fermions. Here e(V) is the differen- 
tial operator such that e(V)e*''''" = ekC^^'^, where ek is the energy spectrum 
of conduction electrons in a specific metal |23]. The attractive interaction 
between electrons is described by a contact pseudo-potential of strength g 
{g > 0). For metals this electron-phonon interaction strength is attractive 
only for conduction electrons near the Fermi surface [TSl [TBI [13 • The chem- 
ical potential /i fixes the number of conduction electrons. 

The Heisenberg equation of motion of the field operator ^/^^(r, t) can be 
immediately derived and reads 

d ^ ... 
^^di^^ = [e(V) - ^]'^f-g ip^^i-^t ■ (2) 

In the BCS theory the interaction term of Eq. ([2]) can be treated within the 
minimal mean- field approximation ip^ipiipf ~ (ipiip^) . In this way Eq. ([2]) 
becomes 

z;i-^t = KV)-/i]^t-A^i , (3) 

where 

A{r,t)=g (^;(r,t)^t(r,t)) (4) 

is the gap function. The condensate wave function of Cooper pairs [Ul [12] 
is instead given by 

S(r,r',t) = (V^^(r,t)^t(r',t)). (5) 

As shown by Yang |l8j, this two-particle wave function is strictly related to 
the largest eigenvalue Nq of two-body density matrix of the system. A^^o gives 
the number of Fermi pairs in the lowest state, i.e. the condensate number of 
Fermi pairs [ISl [13 [IB], and it can be written as 

No = J |S(r,r',t)|2 dhdh' . (6) 

A finite value for the condensate fraction / = Nq/{N/2) in the thermody- 
namic limit N ^ oo implies off- diagonal-long-range-order [22]. 
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3 Gap equation and condensate density 

To investigate the properties of the condensate fraction of electronic pairs we 
adopt the following Bogoliubov representation of the field operator 

^t(r, t) = E (^e^^''--*)6., - :^e-('^--*)6+ ) (7) 

in terms of the anti-commuting quasi-particle Bogoliubov operators fekcr, with 
V the volume of the system and Ek = huk the excitation energies of quasi- 
particles |T5l |16] . 

The thermal averages of quasi-particle Bogoliubov operators are given by 

(^ko-^k'cr') = ^^^^ _^ ^ l^kk'^o-cr' = ?^fc ^kk'^crcr' , (8) 

where /3 = l/{kBT) with the Boltzmann constant, T the absolute tem- 
perature, and fik is the thermal Fermi distribution. 

By using these results, the gap function, Eq. (jlj), becomes 

^ = f E'(l - ^rik)ukVk , (9) 
^ k 

while the condensate number of conduction electrons, Eq. ([6]), satisfies this 
expression [H [19] 

No = - 2n,)\lvl . (10) 

k 

The 'prime' restricts the summation to states within a shell of width huo 
about the Fermi surface. 

To determine the amplitudes Uk and Vk of quasi-particles, one inserts Eq. 
([7j) into Eq. ([3]) and obtains the familiar Bogoliubov-de Gennes equations, 
which give 

where 

^fc = efc - , Ek = + A2 . (12) 
Eqs. ([9]) and (fTOj) can then be written as 

A = ^E'^tanh(^) (13) 
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A^o = E'^tanh^(^) (14) 

where tanh(/3£'fe/2) = 1 — 2nf^. 

In the thermodynamic hmit, where the volume V goes to infinity, J2k can 
be replaced by VJd^k/{2TTf = V f N{^)d^ with iV(0 = Jd^k/{2TTf - 
C,k). In metals the condition huo ^ is always satisfied [16], consequently 
we can use the approximation / N{C,)dC, ~ A^(0) / d^, where 

r d^k 

NiO) = JJ^,5{^^- e,) (15) 

is the density of states at the Fermi surface. In this way the previous equa- 
tions (|T3l) and (|T4l) become 



gN{0) Jo + A' 

1 , , tanh^dve + A") / X 

where tt-q = Nq/V is the density of electrons in the condensate 



3.1 Zero-temperature condensate 



Let us consider first the zero-temperature case (T = 0). From Eqs. ( fT6|) and 
( 1T7|) we get the zero temperature energy gap A(0): 



gNiO) 



In 



D 



+ 



V 



1 + 



D 



A(0)""\ ""A(O) 



(18) 



and the zero-temperature condensate density: 

no(0) = ^iV(0)A(0)arctan(|^: 



(19) 



This expression shows that the condensate density n(0) can be expressed 
in terms of density of states A^(0), energy gap A(0) and Debye energy huD- 
Finally, the zero-temperature condensate fraction /(O) = ?2o(0) /(n/2) is given 
by 



1 A^(0) , , , , huD , 

/ = ^A arctan — 4 

•'^ ^ 2 n ^ ^ ^A ■ 



(20) 
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where n is the density of conduction electrons. 

Under the condition A(0) ^ fiuj^, from Eqs. (llSp we find the famihar 
weak-couphng BCS result 

A(0) = 2hun exp(-^^) (21) 

for the energy-gap order parameter, while the condensate density (fT9l) can 
be written as [211 US] : 

no(0) = ^iV(0)A(0) (22) 

We stress that in many real superconductors the simple BCS theory re- 
ported above is not accurate, and one has to take into account the retarded 
electron-electron interaction via phonons [26J and also the Coulomb repul- 
sion [27] . The results obtained above by using the mean-field BCS theory 
are reliable only in the weak-coupling regime, i.e. for A(0) <C huo, where 
gN{0) < 0.3. Therefore we will continue our analysis of the condensate 
fraction only for a class of superconductors which satisfy this condition. 

Within the free-electrons Sommerfeld approximation, where the energy 
spectrum of conduction electrons has the simple quadratic behavior ek = 
fi^k"^ / (2m*), the free particle density of states A^/ree(0) is related to the 
total density of conduction electrons hj n = 4A^/ree(0)/i/3 and the zero- 
temperature condensate fraction reads /(O) = 37rA(0)/(8/i). To get a better 
estimate, we correct the free electron value of A^(0) by an effective mass 
as obtained from specific heat measurements, i.e. we use the expression 
N{0) = {m*/m)Nfree{0). 

In the first two columns of Tab. 1 we show the zero-temperature con- 
densate density uq and condensate fraction /(O) of simple metals obtained 
from Eqs. ( fT9l) and ( !20l) by using the experimental data of A(0) and ud 
obtained from Ref. 24 (when the comparison is possible, they agree within a 
few percent with those reported in Ref. 25). In the third column we report 
the electron-phonon strength gN{0) calculated with Eq. ( !T8|) knowing /(O). 
The table shows that indeed these simple metals are all in the weak-coupling 
regime. For completeness, in the forth column we insert the theoretical deter- 
mination (see Eq. (12S]1 ) of the critical temperature Tc, which is very reliable 
for these simple metals, when compared with the experimental data. 
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no(0) [10-33 ^-3 




/(O) [10-5] 




gN{0) 


TclK] 


Cd 


4.18 


0.9 


0.179 


0.51 


Zn 


7.72 


1.2 


0.172 


0.79 


Al 


22.4 


2.5 


0.168 


1.15 


Tl 


31.0 


5.9 


0.263 


2.43 


In 


62.1 


10.8 


0.267 


3.46 


Sn 


62.5 


8.4 


0.254 


3.68 



Table 1: BCS predictions for weak-coupling superconducting metals: Hq is 
the zero-temperature condensate density, obtained with Eq. fl22|) : /(O) = 
no(0)/(r2/2) is the zero-temperature condensate fraction; gN{0) is the 
electron-phonon strength, calculated with Eq. f lTSj) . Tc is the critical tem- 
perature from Eq. fl2^ . 

3.2 Finite-temperature condensate 

Let us now investigate the behavior of the condensate density no at finite 
temperature T. Under the condition huo ^ ksTc, which is always satisfied, 
near Tc the energy gap goes to zero according to the power law [Ml [15] 

A(T) = 3.06 ksTc (l - ^) ^ • (23) 

Instead for the condensate density no(T), from Eq. ( IT7I) and the previous 
expression, we find near Tc 

A(TY ( T\ 

noiT) = 0.43 iV(0)-^ = 4.03 N{0)kBTc U - t^) ■ (24) 
kbTc \ TcJ 

For a generic temperature T we solve numerically Eqs. ( fT6l) and (flTl) . The 
theoretical critical temperature T^ obtained from Eq. (JT6l) setting A(Tc) = 0, 
is given by the well-known result [16] 

keTc = 1.13 exp(- ^^^^^^ ) . (25) 

For simple metals the theoretical critical temperature Tc, reported in the last 
column of Table 1, is in good agreement with the experimental one T^^^: the 
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relative difference (^T^^^ — Tc)/T^^p is not large (i.e. within 10%), and for 
some metals (Tl, In, Sn) it is quite small (i.e. within 2%). 

Taking into account Eq. f l20|) . the BCS theory predicts that the zero- 
temperature condensate density in superconductors in the weak-coupling 
regime, can be written as 



This equation resembles the familiar BCS result A(0) = 1.764 /c^Tc for 
the zero-temperature energy gap. 

We stress that the predictions of the BCS theory can be surely improved 
by using the Eliashberg theory [25| [26] . This more sophisticated approach 
will not change the order of magnitude of the numbers in the first two columns 
of Tab. 1 but it could change the last significant figure. 

Coming back to the study of finite-temperature effects, in the upper panel 
of Fig. 1 we plot the condensate density no{T) vs electron-phonon strength 
gN{0) for different values of the temperature T. As expected, by increasing 
the temperature T it is necessary to increase the strength gN{0) to get the 
same condensate density. 

As it happens for the energy gap A(T), one may show that Eqs. ( 1T7|) 
and (fT7|) together also imply that the condensate density may be written as 
its value at T = times a universal function of T/Tc. In the lower panel 
of Fig. 1 we plot the condensate density riQiT) / uq^O) as a function of the 
temperature T/Tc. in the full range of temperatures the numerical results 
(solid line) are reasonably well approximated by (dashed line) 



with a = 3.16 (best fit). 

4 Superconductors vs ultracold atoms 

In metallic superconductors there is an attractive interaction between fermions 
only near the Fermi surface fibi IT6] . On the contrary, as remarked in the in- 
troduction, in superfiuid ultracold two-component Fermi atoms, the effective 
inter-atomic interaction can be made attractive for all atoms of the system 
by using the technique of Fano-Feshbach resonances [121 UHl EHj- This im- 
plies that in the BCS equations for ultracold atoms there is not a natural 



rio(O) = 1.39 N{0)kBT^ . 



(26) 




(27) 
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Figure 1: Upper panel: condensate density nQ{T) vs electron-phonon 
strength gN{0) in a superconductor for different values of the temperature 
T, where N{0) is the density of states and Ed = hun is the Debye energy. 
Lower panel: condensate density no{T) as a function of the temperature T 
in a superconductor with gN{0) = 0.2. Solid line: numerical solution of Eqs. 
(fT6|) and (fT7|) : dashed line: analytical approximation, Eq. (!27|) . 
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ultraviolet cutoff. For attractive ultracold atoms the mean-field BCS theory 
is given by the gap equation and the number equation 



N=}^ivt + 2iui - vi)nu) , (28) 

k 



while the condensate fraction is given by Eq. f lTOj) . But, for ultracold atoms, 
in these equations the sum over momenta is no more restricted within a thin 
shell around the Fermi surface. As well known, due to the choice of a contact 
potential, the gap equation ^ diverges in the ultraviolet. This divergence 
is logarithmic in two dimensions (2D) and linear in three dimensions (3D). 

In 3D, a suitable regularization is obtained by introducing the inter- 
atomic scattering length ap via the equation 

m 1 1 1 

'^l^P^~~g^v\2r,' ^^^^ 

where = h'^k'^/{2m) with m the atomic mass, and then subtracting this 
equation from the gap equation O [6l [7]. In this way one obtains the 3D 
regularized gap equation 



In 2D, quite generally the bound-state energy exists for any value of the 
interaction strength g between two atoms [29]. For the contact potential 
the bound-state equation is 

and subtracting this equation from the gap equation [10], |29] one obtains a 
2D regularized gap equation 




tanh(/3Efc/2) 
2Ek 



. (32) 



The number equation fl28|) and the renormalized gap equation (l30l) (or Eq. 
fl32|) in 2D) are the so-called generalized BCS equations, from which one 
determines, for a fixed value of the temperature T and the average number 
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of atoms A^, the chemical potential /u(T) and the gap energy A(T) as a 
function of the scattering length ap (or of the bound-state energy in 2D). 
The extended BCS equations can be applied in the full crossover from weak 
coupling to strong-coupling [151 l2H]- In 3D, the crossover is from the BCS 
state of weakly-interacting Cooper pairs (with l/ap ^ —1) to the Bose- 
Einstein Condensate (BEC) of molecular dimers (with l/ap ^ 1) across the 
unitarity limit [l/ap = 0) [6J. In 2D, there is a similar BCS-BEC crossover 
by increasing the value es of the bound-state energy [101 II]- 

At zero-temperature, by using the continuum limit X^k ~^ ^ I (27r)^ / (i^k — )• 
V j (27r^) / fc^cifc, the 3D condensate density ( ITOj) has a simple analytical ex- 
pression [I]. The 3D density of states is N{i) = {2m/h'^f/^y^lTJl/{47r^) 
and the 3D condensate density is given by 



m 



3/2 



MO) = ^A(0) 



3/2 



l + (33) 



87ih^ \ A(0) \l A(0) 

In the 3D BCS regime {l/ap ^ —1), where /i(0)/A(0) ^ 1 and the size 
of weaklybound Cooper pairs exceeds the typical interparticle spacing kp^, 
/i(0) approaches the non-interacting Fermi energy ep = h^k'^/ {2m) with 
kp = (37r^n)^/'^ and there is an exponentially small energy gap A(0) = 
Se~'^ep eyi\i {n / {2kpap)). In this weak-coupling regime the 3D condensate 
density becomes [1] 

"o(0)4iV(0)A(0)^|L„e.p(5jl-). (34) 

Notice that this is formula is similar to Eq. f l2^ of weak-coupling supercon- 
ductors (here ai? < 0), but the behavior of A(0) is quite different. 

In 2D, the density of states is constant and reads N{^) = N{0) = 
(2m/^^)/(47r). The zero-temperature 2D condensate density is easily obtained[l] 

as 

no(0) = ^Ar(0)A(0) (^^ + arctan {^)^ , (35) 
while the zero-temperature 2D energy gap is given by the implicit formula 



(36) 
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From these equations, in the 2D BCS regime (0 < ^ ^f) where /i(0)/A(0) ^ 
1 one finds exactly Eq. ( 122]) . but here the energy gap A(0) depends on the 
Fermi energy ep and the bound-state energy according to the formula [29] 

A(0) = , (37) 

while the chemical potential is fi{0) = ep — £5/2. It is not surprising that in 
the BCS regime the condensate density of 2D superfluid atoms is formally 
equivalent to the Eq. fl22l) we have found for weak-coupling superconductors. 
In fact, to obtain Eq. fl221) we have used the approximation / N{C,)dC, — 
N{0) J dC, that is exact in the strictly 2D case, and the condition A(0) <^ huu 
which implies that the upper limit of integration is practically +00. 

In the previous section we have shown that the BCS equations can be 
used to determine the (quite small) condensate fraction of superconductors 
only in the weak-coupling regime. Instead, the extended BCS equations 
have been used in recent papers [H El |2l S] to get the condensate fraction of 
ultracold atoms in the full BCS-BEC crossover. The theory predicts that in 
the crossover the zero-temperature condensate fraction grows from zero to 
one. Two experiments [HI [T3] have confirmed these predictions for the 3D 
superfiuid two-component Fermi gas. 

5 Conclusions 

In this paper we have studied, within the mean-field BCS theory of super- 
conductors, the condensate of electronic Cooper pairs at zero and finite tem- 
perature showing the crucial role played by the Debye frequency and by the 
electron-phonon interaction. We have found that the zero-temperature con- 
densate fraction /(O) of weak-coupling metals is quite small (~ 10~^) and 
the condensate density increases in metals with higher critical temperature 
Tc, according to the law /(O) = 1.39 N{0)kBTc, where iV(0) is the density 
of states at the Fermi energy. As discussed by Chakravarty and Kee [2T] . 
the spin-spin correlation function depends significantly on the condensate 
density and magnetic neutron scattering can provide a direct measurement 
of the condensate fraction of a superconductor. In the next future our BCS 
predictions, which are meaningful for weak-coupling superconductors, could 
be experimentally tested. In the last part of the paper we have shown sim- 
ilarities and differences between metallic superconductors and atomic Fermi 
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vapors in the determination of the condensate fraction by using the mean- 
field BCS theory and its extension in the BCS-BEC crossover. 

The author thanks A.J. Leggett and F. Toigo for useful comments and 
critical remarks. 
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